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Abstract 

A two-step doping process, magnetic followed by charge or vice versa, is required to produce 
insulating massive surface states in topological insulators for many physics and device applications. 
Using first-principles calculations, we demonstrate here simultaneous magnetic and hole doping 
achieved with a single dopant, carbon, in Bi2Se3. Carbon substitution for Se (Cse) results in an 
opening of a sizable surface Dirac gap (53-85 meV), while the Fermi level (Ep) remains inside the 
bulk gap and close to the Dirac point at moderate doping concentrations. The strong localization 
of 2p states of Cs e favors spontaneous spin polarization via a p-p interaction and formation of 
ordered magnetic moments mediated by the surface state. Meanwhile, holes are introduced into 
the system by Cs e - This dual function of carbon doping suggests a simple way to realize insulating 
massive topological surface states. 

PACS numbers: 68.35.Dv, 73.20.At, 71.15.Mb 
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The second generation 3D topological insulators (TIs), Bi 2 Se3 and Bi 2 Te3, are a class of 
time-reversal-invariant materials characterized by an insulating bulk state and conducting 
surface state consisting of a single Dirac cone at the T point (see Refs. [l, 2] and references 
therein). Such massless surface states are protected by the time-reversal symmetry (TRS) 
and are immune to surface disorders such as defects Shi]]. Many striking physics and 
device applications of TIs have been proposed, such as quantum anomalous Hall effect 
, magnetic monopole imaging |8||, topological contribution to the Faraday and Kerr 
effects jjj, inverse spin-galvanic effect 10], optical injection of spins 11] and Tl-based p-n 



junctions 



12 



13] . To enable the above applications, it is necessary to open a surface energy 



gap as well as keep the Fermi energy inside the bulk gap|9|. It was proposed that doping 
TIs by magnetic transition metals (TM), such as Fe, Mn, Cr and Co, could break the TRS 
and open a surface gap|3|, which were commonly used experimentally now 

aai 

14H22J. In 



the experiment, the existence of anionic vacancies in the as-growth Bi 2 Se 3 or Bi 2 Te3 thin 



film lead to a typical n-type material |15|, |23j, |24j. To move the Fermi level inside the 
bulk gap from the bulk conduction band, the excess electrons must be compensated (vi a 



y 



hole doping) which can be achieved by doping with divalent cations (Mg 2 " 1 " or Ca 2+ ) 14], 



or chemical molecules with strong e 
such as ternary Bi .o8Sbi.g 2 Te3 {Jl^ 



ectron accepting ability (N0 2 , CO) |l5l . |23I| . Alloying, 



24] and quaternary Bii. 5 Sbo.5Tei.8Sei. 2 [25], is another 



possible approach to compensate the excess electrons but it is difficult to accurately control 
the compositions of these compounds [3, S, 25]. Nevertheless, it remains a challenge to find 



a reliable yet simple scenario to produce out-of-plane magnetization, a stable long-range 
magnetic order and proximity of Fermi level to the Dirac point. These 3 conditions are 
necessary in order to experimentally observe the Dira c gap and quantum anomalous Hall 
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-|30|. 



effect using the approach of magnetic doping {3], 7, 

The spontaneous spin polarization and local moments in carbon-doped ZnO, a diluted 
magnetic semiconductor (DMS), has been successfully demonstrated both experimentally 



and theoretically 
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32] . Since O and Se have similar electronic configurations and similar 



chemical properties, and C has a quite similar electronegativity (2.55) but two more valance 
holes compared to Se, one can expect that substitution of C for Se would compensate both 
the Se vacancies and excess electrons in the as-growth Bi 2 Se3. This motivated us to consider 
carbon doping as a mean to achieve massive topological surface states (TSS) in Bi 2 Se3. 
In this Letter, we report results of our investigation on the topological surface state 



of carbon doped Bi 2 Se3 using first-principles electronic structure calculations. It is found 
that substitution of C for Se in Bi2Se3 introduces local moments as well as holes. Carbon 
doping can lead to simultaneous opening of the Dirac gap up to 85 meV and pinning of the 
Fermi level inside the bulk energy gap. It is striking that magnetic doping and ordering can 
be achieved using a nonmagnetic dopant (C). We discuss the origin of magnetic moment, 
mechanism of magnetic coupling, and cause of magnetocrystalline anisotropy in carbon 
doped TIs system, and compare them with TM-doped TIs and C-doped DMSs. The effects 
of carrier and impurity concentration on the magnetic state and Dirac gap are also discussed. 



First-princip 



great success 



es calculations have been widely used to predict topological insulators with 



331 ] . Our First-principles calculations were carried out using the VASP code 34 1 



with the projector-augmented-wave potentials and the Perdew-Burke-Ernzerhof generalized 



gradient approximation 35[ for electron exchange-correlation functional. The experimental 
lattice constants of Bi2Se3 were used in our calculations [ssj] . The Bi2Se3 thin film was 
modeled by a slab of six quintuple Se-terminated layers (QLs). A vacuum layer of more 
ihan 30 A thick was included to ensure a negligible interaction between neighboring slabs 
23), 0]- AH structures were relaxed until the force on each atom is smaller than 0.01 eV/A . 
The cutoff energy was taken to be 400 eV and /c-point meshes of 7x7x1, 5x5x1 and 
3x3x1 were used for the 2 x 2, 3 x 3, and 5x5 surface unit cell, respectively. Different 
surface unit cells were used to study effect of doping concentration on the magnetic and 
electronic properties. The SOI was included in the topological surface state calculations. 
Since Se vacancies are easily formed in Bi 2 Se3 (Bi-rich growth condition experimentally) 



{14I l23l l38l l39j . we first considered a Se vacancy (Vs e ) on the Bi2Se3 surface which was 
modeled by removing a Se atom from the top quintuple layer in the supercell. Due to the 
layered structure of Bi 2 Se3, there are three inequivalent Vs c sites, with Vs c at the top (1), 
middle (2) and bottom (3) of the Bi 2 Se3 layer, as shown in Fig. la. The calculated vacancy 
formation energies are 0.661 eV, 0.475 eV and 0.158 eV for sites 1-3, respectively. The 
vacancy formation energy is defined as E = E^ef — E pr[n + n/ise, where E^d and E pTin are 
total energies of the film with and without the Se vacancy, fi^e is the chemical potential 
of atomic Se, and n is number of Se vacancies in the supercell. In order to simulate the 
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experimental Bi-rich growth condition |l_4j, |23|, |38|, |39], the chemical potential of atomic Bi 



(/iBi) was taken to be the total energy per atom in Bi bulk phase, while [i$ c is obtained 
from 2// B i + 3/ise = E Bi2Se3 , with E-Q i2 a, e . A being the total energy of bulk Bi 2 Se 3 per molecular 



formula. The calculated formation energies indicate that Se vacancy is more likely to form at 
site 3 in Fig. la. The band structure of the 6QL Bi2Se3 thin film with a single Se vacancy at 
site 3 is shown in Fig. lb. As can be seen, the topological surface state with the Dirac-cone 
characteristics is preserved within the bulk band gap in the presence of the Se vacancy, but 
the Fermi level Ep is located in the bulk conduction band (BCB), almost 0.5 eV above the 
Dirac point, due to the anionic Se vacancy. This is in good agreement with experimental 



observations (inset of Fig. lb) [39 1 



We next consider filling of the Se vacancy by carbon. Theoretically, it should not be 
difficult to fill the Se vacancy in Bi 2 Se3 with carbon. This is because carbon has the same 
electronegativity (2.55) as Se, but much smaller atomic radius (70 pm). Furthermore, the 
large interstitial space in the Bi 2 Se3 crystal lattice, due to the large covalent radius of Bi (146) 
and Se (116), is able to facilitate diffusion of carbon to Vs e sites on the Bi2Se3 surface. As the 
doping concentration increases, carbon can further replace selenium during high temperature 
annealing [l^]. Figure lc shows the band structure of Bi 2 Se3 in the same configuration as 
in Fig. lb but the Vs e is filled by a carbon, resulting in a carbon substitutional doping. It 
shows that carbon doping leads to the opening of a sizable surface Dirac gap (53 meV), 
which is comparable to that in Fe or Mn doped Bi 2 Se3 [l^, [lsj]. Such a large Dirac gap is 
desirable for room temperature device applications. In addition, the Fermi level is drawn 
closer to the surface gap from the bulk conduction band due to hole doping introduced by 
Cg~ 2 - (see Fig. la). This indicates that we can realize insulating massive TSS in one step, 
in contrast to the two-step process by TM and carrier doping. Experimentally, a one-step 
doping procedure would be much preferred over the two-step process because it means easy 
control of experimental environment and parameters. 

In order to understand how a nonmagnetic element such as carbon can open a Dirac 
gap and tune the Fermi level simultaneously, we calculated the projected density of states 
(DOS) of carbon doped Bi 2 Se 3 thin film and present the results in Fig. 2. A strong coupling 
between the 2p states of the dopant (C) and the 4p/Qp states of the host (Se/Bi) is clearly 
seen at the valence band maximum (VBM) and conduction band maximum (CBM). This 
p-p interaction is essentially a result of quantum-mechanical level repulsion, which splits 
the majority and minority 2p states of carbon and "pushes" the minority 2p states up 
inside the bulk band gap. The exchange splitting energy between the majority and minority 
spin t 2 states at the zone center [Aer = e(t|) — e(t 2 )] is substantial (~0.80 eV). Further 



examination on the projected p states of carbon [see inset of Fig. 2] shows the p x and p y 
orbitals are degenerate. The spin-down p XyV states are partially occupied, while the spin- 
down p z states are completely empty. The energy difference (Ae n ) between pj. y and p\. 
is around 0.66 eV. Such a splitting is due to the change of carbon bonding environment, 
especially in the z direction, under the tetrahedral crystal field. The symmetry and wave 
function of the impurity state (p-like £2) are similar to those of the VBM and CBM of Bi2Se3 
which consists mainly of anion Se 4p and cation Bi 6p orbitals (without SOI) [see Fig. 2] 
or cation 6p and anion 4p orbitals (with SOI) due to band inversion. Therefore, a strong 
p-p coupling between the impurity state and valence-band state is allowed near the Fermi 
level, resulting in the spin polarization and formation of local moments. Such local moments 
are quite stable and the spontaneous spin-polarization energy is 175 meV. Our calculations 
also show (see discussion on impurity concentration effect below) that these local moments 
favor an out-of-plane long-range ferromagnetic coupling on the surface of Bi 2 Se3, mediated 
by the surface state. The carbon valence orbitals favor more stable fourfold (T = a\ + 1 2 ) 
splitting under the tetrahedral crystal field, forming C 4 ~ (see Fig. 1), which provides two 
more holes compared to Se 2- . Therefore, Cs e can introduce holes as well as local moments. 
This duality of carbon doping leads to the simultaneous topological surface state and Fermi 
level tuning in this intrinsic nonmagnetic topological system. In topological insulators, light 



elements are often used to tune the surface hole doping by modifying their concentration [40j. 
Experimentally, the carrier and impurity density are two key parameters. We therefore 
discuss their effects on the stability and Fermi level tuning in carbon-doped Bi2Se3. 

Carrier effect — Carriers play critical roles in magnetically doped TIs, in both bulk and 
surface. They mediate the long-range magnetic order 18| and control the Fermi level as well 
. Electrons in the system also affect the localization of magnetic moments, especially 
in systems with delocalized p band and hole-induced magnetic system. Here, we consider 
only the effect of electrons on the local moments and the Dirac gap in C-doped Bi 2 Se3. As 
shown in Fig. 3, C substitutional doping at the anion site induces a magnetic moment of 
~2/xg per anion. If excess electrons are introduced into the system which is possible in an 
experimental environment, they would quench the magnetic moments, the Dirac gap can 
be closed if the carrier concentration is sufficiently high. We wish to point out, however, in 
our calculation the additional electrons are assumed to be homogenously distributed in the 
whole system, rather than localized at or around the carbon dopants. 



Impurity concentration effect — The impurity concentration can affect both the degree of 
localization of magnetic moments and the hole density because the carbon dopants introduce 
local moments as well as holes into the system. Figure 4 shows the dependence of the local 
moment and Dirac gap on the impurity concentration which is given in terms of distance 
between impurities, by assuming a uniform dopant distribution. It is noted that the inter- 
impurity distance in the figure is discontinuous. A high impurity concentration enhances 
the magnetic coupling and introduces more holes, but too large an impurity concentration 
would quench the local moment and reduce the size of the Dirac gap. This can be under- 



stood based on the phenomenological band-coupling model [41]. Under a very low impurity 
concentration, the impurity states are localized and isolated, which is insufficient to open a 
Dirac gap. At a moderate doping concentration, the 2p states with the same spin couple to 
each other but remains localized (FM-I phase), and a Dirac gap opens. The localized im- 
purity bands become broader and delocalized (FM-II phase) as the impurity concentration 
increases further. This eventually leads to a charge transfer from the majority spin state to 
the minority spin state, and a reduced local moment and exchange splitting (Fig. 4), and a 
corresponding reduction of the Dirac gap. If the impurity concentration is sufficiently high, 
the much reduced exchange splitting would lead to an AFM ground state (AFM phase), 
even though the energy gain of the FM phase is usually larger than that of the AFM phase 
in C-doped TI due to the second-order nature of the superexchange interaction of the AFM 
phase. Therefore, a moderate carbon doping concentration is necessary to open a sizable 
Dirac gap and to tune the Fermi level close to the Dirac point. 

C-doped TIs vs. TM-doped TIs — There are several important differences between the 
C-doped and Fe/Mn-doped TIs due to the distinctly different orbitals (2p vs. 3d) of the 
dopants. First, the p orbitals of C are usually fully occupied in ionic states, leaving no room 
for unpaired electrons compared to the d orbitals of transition metals. Second, the spin-orbit 
and hyperfine interactions in C are considerably weaker compared to that in TMs since they 
scale as the fourth power of the atomic number. Consequently, it is possible to preserve 
spin-coherence over time and distance in C-doped TIs much longer than in TM-doped TIs, 
especially for the massive particles in the presence of a Dirac gap. Finally, concerning the 
mechanism of magnetic ordering in TIs, Liu et al., proposed mediation of the long-range 



ferromagnetism by surface states through the RKKY interaction in TM-doped TIs 



26. 



while Chang et al. suggested the p-d van Vleck mechanism (localized p valence electrons) 



for the long-range magnetic order in TM-doped TIs [22|. Clearly, the ferromagnetism in 
carbon doped Bi2Se3 challenges our understanding of magnetic ordering in such systems 
because there are no d states near either conductance or valence bands in C-doped TI 
system (see Fig. 2). The above s-d RKKY or p-d van Vleck mechanism cannot be applied 
here. The ferromagnetic double-exchange mechanism can produce large spin moments, but 
it is a short-range interaction that requires mixed valence, i.e., 2p n 2p n+1 . However, 
there is no evidence that mixed valence occurs in C-doped Bi2Se3 based on our magnetic 
moment calculations. We propose that the spontaneous spin polarization is induced by the 
p-p interaction and the magnetic exchange coupling in such a system is mediated by the 
surface states 22|. Further investigation is required to clarify the magnetic exchange in 
C-doped Bi 2 Se3 as well as TIs doped by other 2p light elements in general. 

C-doped TIs vs. C-doped DMSs — We also wish to point out an important difference 
between magnetism in C-doped TIs and C-doped DMSs, which is essentially due to the 
different atomic orbitals in the host materials (3d in DMS vs. 4/ in TIs). Bismuth surface 



shows a strong SOI due to its high atomic numbers 42J. A large orbital moment is observed 



on the surface 



421 ] . which is driven by a localized spin through SOI. Such an orbital moment 



gives rise to an effective field that is responsible for the large surface magnetic anisotropy in 
C-doped Bi 2 Se 3 (around 10 meV per C), which is comparable to that in Fe-doped Bi 2 Se 3 [3] 



and Co-doped Bi 2 Se3 (43] . The magnetic order in C-doped TIs can be expected to survive 
up to room temperature since thermal energy is insufficient to destroy such a magnetic 
anisotropy. In contrast, the SOI in DMSs is much weaker since the host materials are 
typically GaAs, ZnO et al. Therefore, the Curie temperature of most DMSs is below room 



temperature despite a few exceptions such as C-doped ZnO 



311 ] . C-doped TIs and C-doped 



DMSs are also different in their dependence on carriers. The insulating FM phase of C-doped 



TIs can be carrier-independent for the QAH effect and spin Hall device applications [22 ]. 
This is different in C-doped DMSs where carriers (hole) of certain density are indispensable 
for the long-range FM coupling between local moments and for spin field-effect transistor 



(FET) applications 3jJ 



In summary, we propose carbon substitutional doping at Se site as a one-step approach 
to achieve massive topological surface states in Bi 2 Se3, which is more effective compared 
to the usual approach of magnetic impurity doping followed by chemical charge doping 
or vice versa. This is possible because carbon doping simultaneously introduces localized 



spin moments and holes. Such spin moments are long-range ordering by p-p interaction. 
Meanwhile, holes introduced by carbon doping pull the Fermi level inside the bulk gap 
from the bulk conduction band. We expect this simple approach of producing massive 
topological surface states to promote experimental studies in fabricating C-doped TIs or 
similar materials, and eventually lead to simpler experimental procedure and robust TIs for 
broader device applications, such as CMOS-like, spin logic and next-generation low-power 
devices utilizing the charge on/off states (because of opened gap), spin up/down states 
(because of topological surface) and Hall conductance (because of anomalous Hall effect). 
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a b c 




sp orbitals splitting in a 

tetrahedral crystal field ( V Se ) ( C Se ) 



FIG. 1. (Color online) (a) Different anionic vacancy /doping sites on the surface quintuple layer of 
a Bi2Se3 slab, and schematic diagram of carbon substitutional doping at site 3 and splitting of its 
orbital under the tetrahedral crystal field, (b) Calculated band structure of 6QL Bi2Se3 thin films 
with a Se vacancy defect. The inset at the upper-right corner is the experimental band structure 
of 6QL Bi2Se3 thin films [ref. 39]. (c) Calculated band structure of 6QL Bi2Se3 thin films with 
one carbon substitutional doping. The projected pristine bulk band structure to the 2D Brillouin 
zone are shown in the background (blue). The surface topological states are highlighted by the 
bold purple lines. The insets (lower-right corners) are schematic illustrations of band structures of 
Bi2Sc3 with Se vacancies (as-growth) and substitutional carbon doping at Se site (Cs e ), respectively. 
The Dirac energy level is labelled at the Dirac point without a gap and the bottom of the surface 
"conduction" band with a gap, respectively. 
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FIG. 2. (Color online) DOS (without SOI) of carbon doping in the Bi2Se3 thin film. The up panel 
is majority spin and the down panel is minority spin. The inset is the projected DOS on carbon 
2p orbitals. 
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FIG. 3. (Color online) Dependency of spin moment and Dirac gap on the carrier (electron) con- 
centration. 



13 



concentration 



2.0 



m 


1.5 






-t— ' 

c 




CD 




E 


1.0 


o 




E 




c 




Q. 
CO 


0.5 








o 



coupled I — i 
localized [ — I 



FM-I FM-II AFM 
~ J coupled 
\ localized 



, isolated 
-| I lo calized 




coupled \\ 
del ocalized Q- 



20.7 8.3 7.2 4.1 3.3 

o 

Impurity distance (A) 



1*1 

U t i 

Isolated 

r> 

T 4 T i 
Isolated 



FM-I 



FM- 



AFM 



FIG. 4. (Color online). (a) Dependence of magnetic moment and Dirac gap on the inter-impurity 
distance, which is used to indicate impurity concentration, assuming a uniform impurity distribu- 
tion, (b) Schematic illustration of the possible magnetic coupling in carbon doped TIs. The left 
figure shows the DOS for two isolated carbon dopants, and the right figure shows the coupling 
when the two carbon dopants are in close proximity at a higher doping concentration. The dashed 
lines denote the Fermi level. 
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